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Abstract. In this article we present some intcgrability conditions for partial difference 
equations obtained using the formal symmetries approach. We apply them to find in- 
tegrable partial difference equations contained in a class of equations obtained by the 
multiple scale analysis of the general multilinear dispersive difference equation defined 
on the square. 



1. Introduction 

DL met Marcos Moshinsky for the first time on his arrival in Mexico in February 1973. 
He went there with a two years fellowship of the CONACYT (Consejo Nacional de Ciencia 
y Tecnologia), the Mexican Research Council. DL visit was part of his duties as military 
service in Italy and was following a visit to Mexico in the summer of 1972 by Francesco 
Calogero, with whom he graduated at the University of Rome La Sapienza in the spring of 
1972 with a thesis on "Computation of bound state energy for nuclei and nuclear matter 
with One-Boson-Exchange Potentials" [5]- 

The bureaucratic process, preliminary to this visit, was long and tiring and it was suc- 
cessful only for the strong concern of DL uncle, Enzo Levi, professor of Hydrulics at UN AM 
[S] and for the help of Marcos Moshinsky who signed a contract for him which, as DL uncle 
wrote, was ready to substain with his own personal funds if no other source could be found. 

During DL stay in Mexico he was able, with the help of Marcos, to appreciate the plea- 
sures of research and to start his career in Mathematical Physics. With Marcos Moshinsky 
he published, during the two years stay at the Instituto de Fisica of UN AM, 4 articles 
[3 HU [7[ [5] but he also collaborated with some of the other physicists of the Institute 
[HEHE2]. This visit was fundamental to shape DL future life, both from the personal 
point of view and from the point of view of his career. 

After 1975 DL met Marcos Moshinsky frequently both abroad and in Mexico. When 
DL was in Mexico, Marcos Moshinsky always invited him to his home to partecipate to his 
family dinners. As a memento of Marcos, let us add a picture of his with Pavel Winternitz, 
Petra Seligman and Decio Levi, taken in September 1999 in front of the Guest House of 
the CIC-AC where he stayed overnight after a seminar given at the Instituto de Fisica of 
UNAM, Cuernavaca branch. 

Symmetries have played always an important role in physics and in the research of 
Marcos Moshinsky. This presentation shows how crucial can be the notion of symmetry in 
uncovering integrable structures in nonlinear partial difference equations. 

The discovery of new integrable Partial Difference Equations (PAE) is always a very 
challenging problem as, by proper continuous limits, we can obtain integrable Differential 
Difference Equations (DAE) and Partial Differential Equations (PDE). 

A very successful way to uncover integrable PDE has been the formal symmetry approach 
due to Shabat and his school in Ufa [32] . These results have been later extended to the case 
of DAE by Yamilov |31j . a former student of Shabat. 

Here we present some results on the application of the formal symmetry technique to 
PAE. 

The basic theory for obtaining symmetries of differential equations has been introduced 
by Sophus Lie at the end of the nineteen century and can be found together with its 
extension to generalized symmetries introduced by Emma Noether, for example in the book 




by Olver [55] . The extension of the classical theory to PAE can be found in the work of 
Levi and Winternitz [H i2 HH] . 

Here in the following we outline the results of the geometric classification of PAE given 
by Adler, Bobenko and Suris |2j with some critical comments at the end. Then in Section 2 
we derive the lowest integrability conditions starting from the request that the PAE admit 
generalized symmetries of sufficiently high order. In Section 3 we show how the test can be 
applied. Then we apply it to a simple class of PAE obtained by the multiple scale analysis 
of a generic multilinear dispersive equation defined on the square. 

1.1. Classification of linear affine discrete equations. Adler, Bobenko and Suris [5] 
considered the following class of autonomous PAE: 

Ut+i.j+i = F(u i+1 j,Uij,Uij+i), (1) 

where i,j are arbitrary integers. Eq. ([I]) is a discrete analogue of the hyperbolic equations 

u xy = F (u x ,u, u y ). (2) 

which are very important in many fields of physics. Up to now the general equation ^ has 
not been classified. Only the two particular cases: 

u xy =F(u); u x = F(u,v) , v y = G(u,v), 

which are essentially easier, have been classified by the formal symmetry approach |321 133j . 

The ABS integrable lattice equations are defined as those autonomous affine linear (i.e. 
polynomial of degree one in each argument, i.e. multilinear) partial difference equations of 
the form 

£(u , ,ui fi , u 0< i, tti.i; a, f3) = 0, (3) 

where a and j3 are two constant parameters , whose integrability is based on the consis- 
tency around a cube (or 3D-consistency) . Here and in the following, as the equations are 
autonomous, and thus translational invariant, we skip the indices i and j and write the 
equations around the origin. 

The main idea of the consistency method is the following: 
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Figure 2. A square lattice 
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Figure 3. Three-dimensional consistency 



(1) One starts from a square lattice and defines the three variables Uij on the vertices 
(see Figure 2). By solving £ — one obtains a rational expression for the fourth 
one. 

(2) One adjoins a third direction, say k, and imagines the map giving "1,1,1 as being 
the composition of maps on the various planes. There exist three different ways 
to obtain "1,1,1 and the consistency constraint is that they all lead to the same 
result. 

(3) Two further constraints have been introduced by Adler, Bobenko and Suris: 

• D4- symmetry: 

£("o,o,"i,o,"o,i,"i,i;a,/3) = ±£("0,0, "0,1, "1,0, "1,1; P, ex) 

= ±£("i,o,"o,o,"i,i,"o,i;a,^)- 

• Tetrahedron property: "1,1,1 is independent of "0,0, o- 

(4) The equations are classified according to the following equivalence group: 

• A Mobius transformation. 

• Simultaneous point change of all variables. 

As a result of this procedure all equations possess a symmetric (in the exchange of the 
first to the second index) Lax pair, Backlund transformations etc. . Thus the compatible 
equations are, for all purposes, completely intcgrable equations. 
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The ABS list read: 

(HI) (u ,o - (ui,o - «o,i) - " + /3 = 0, (4) 

The potential discrete KdV equation [T3 [53] 

(H2) (u ,o - ui,i)(«i,o - uo,i) + (P ~ a )( u o,o + "i.o + uq,i + ux,x) ~ 

-a 2 + /3 2 = 0, 

(H3) a(uo,oWi,o + uo,xux,x) - /8(ito,oWo,x + tti,o«x,x) + ~ £ 2 ) = °> 

(5) 

(Ql) a(uo,o - uo,i)(«x,o - wx.x) _ P(uo,o - «i,o)(«o,x - "14) + 

+5 2 a/3(a — /3) = 0, The Schwarzian discrete KdV equation [TBI |2"5] 

(Q2) a(u ,o - «o,i)(«i.o - "x,x) - P{ u o,o - «x,o)(«o,x - «x,x) + 

— /8)(uo,o + «x,o + "o,x + u i,i) — Of/8 (a — /3)(a 2 — a/? + /3 2 ) = 0, 
(Q3) (/3 2 - a 2 )(u ,o'"x,x + «x,o«o,i) + /3(a 2 - l)(uo,o«x,o + wo,x«x,x) - 

-a(/? 2 - l)(uo, Uo,i + «i,o«i,i) , -a — 1 = 0, 

4ap 

(Q4) a tto,oWi,otto,xMx,x + 

+ai(uo,oWx,o'"o,x + «i, o"o,xUx,x + "0,1^1,1^0,0 + ui,i«o,o*ti,o) + 

+a 2 (M ,0Ml,l + WX,0«0,l) + O2(«0,0«X,0 + "O.X^X.l) + 

+a 2 (u , 0^0,1 + Wx,0«X,l) + a 3("0,0 + «X,0 + "0.1 + "X,l) + «4 = 0, 

where the seven parameters Oj's in (Q4) are related by 3 equations. 

By a proper limiting procedure all equations of the ABS list are contained in eq. (Q4) [2~4"| . 
The symmetries for the discrete equations of the ABS list have been constructed [28, 27 and 
are given by DAE, subcases of Yamilov's discretization of the Krichever-Novikov equation 
(YdKN) HZH3T]: 



duo R(ui,uo,U-i) 
de ui — U—i 

where 



R(ui,u , u-i) = AqUiU-x + B {ui + u~i) + C , 



A = ci«o + 2c 2 u + c 3 , 



Bq = c 2 ul + c 4 u + c 5 , 



Co = c 3 u + 2c 5 u Q + c 6 . 

It is immediate to see that by defining = Uij and Vi = Uij + i, the equations of the ABS 
list are nothing else but Backlund transformations for particular subcases of the YdKN 
[T3l 117] . The ABS equations do not exhaust all the possible Backlund transformations for 
the YdKN equation as the whole parameter space is not covered [JTIISQ]. Moreover, in the 
list of integrable DAE of Volterra type [3JJ, there are equations different from the YdKN 
which may also have Backlund transformations of the form ([I]). So we have space for new 
integrable PAE which we will search by using the formal symmetry approach. An extension 
of the 3D consistency approach has been proposed by the same authors [3] allowing different 
equations in the different faces of the cube. However in this way ABS were able to provide 
only examples of new integrable PAE but not to present a complete classification scheme. 

2. Construction of Integrability Conditions 

We consider the class of autonomous PAE 

ui,i = /o,o = -F(wi,o,wo,o,Wo,x) (d Ul . a F, d Uoo F, d U01 F) ^ 0. (6) 

Introducing the two shifts operators, T\ and T 2 such that TiUij — itj+i,^, T 2 Ui j — mj+i, 
it follows that the functions related among themselves by eq. (rol and its shifted 

values 

u i+i,j+i = TiT^fofi = fij = F(u i+ ij,Uij,Uij + i). 
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So, the functions Uij are not all independent. However we can introduce a set of independent 
functions Uij in term of which all the others are expressed. A possible choice is given by 
(u.i Q, u j), for any arbitrary i,j integers. 

A generalized symmetry, written in evolutionary form, is given by 



dt 



uo.o = .9o,o = G(u n ,o, Un-i,o, ■ ■ ■ ,"n',o,uo,fc,wo,fc-i, ■ ■ ■ ,Uo,k'), n > n' , k > k'. (7) 



where t is the group parameter. By shifting, we can write it in any point of the plane 

^ U i,j = ^1^2 #0,0 — 9i,j — G(u i+n j, . . . , Ui+n'j, Uij+ki ■ ■ ■ , u i,j+k')- 

In term of the functions gij we can write down the symmetry invariant condition 

d/o,o" 



9i,i 



dt 



0. 



(8) 



«i,i=/o,o 



i.e. (71,1 = (gifldu 10 +go,od UOQ + go,i9 UQ J/o.o- This equation involves the independent 
variables (iti.o, u j) appearing in g ,o shifted to points laying on lines neighboring the axis, 
i.e. (ui.i, ui.j). For those function we can state the following Proposition [50], necessary to 
prove the subsequent Theorems: 

Proposition 1. The functions Ui^U\,j have the following structure: 

i > : u it i = ^,1(^,0,^-1,0, • • • ,"i,o, "0,0,^0, 1), ^u,, M i,i = r i _1 /«i, ; 



d Ui0 Ui,i = ~T{ 



f , 



i < : Ui,i = Ui^Uifi, u i+ ifl, . . . ,u_i,o, wo,o, "0,1), 
j >0: uij = uij(v, li o,u fl,uo il ,...,Uo,j-i,uo,j), d Uoj u ld = T^~ L f u - 
j <0: uij = ui,j(ui, ,uo,o,"o -1, ■ ■ ■ , woj+i, wo,j), 'I- . " 1 ,j =-T| 7 ^ 1 . 



(9) 



In eq. |9| and in the following, f UiJ = and # Ujj 



dgo.o 



If a generalized sym- 



metry of cnaracteristic function 50,0 depends on at least one variable of the form Ui t o, 
then (g Un0 ,9u n , ) 7^ 0, an< i the numbers n,n' are called the orders of the symme- 
try. The same can be said about the variables uqj and the corresponding numbers fc, k' if 

{9u , k ,9u 0th ,) ^ 0. 

Now we can state the following Theorem, whose proof can be found in [20] : 

Theorem 1. If the PAE u± t i = F possesses a generalized symmetry then the following 
relations must take place: 



n > 0, 
ri < 0, 
k > 0, 
k' < 0, 



(T? - l)log/ Ul . - (1 - T 2 )Tilog 5u „, , 



(Tf -l)log 



fuo,i 



(1-T 2 ) log 



(If - 1) log / Uo>1 = (1 - T X )T 2 log 5uo _ fc , 
(T 2 fc '-l)log^ = (l-ri)log ffuofc; . 



(10) 

(11) 
(12) 
(13) 



As 



TJ n -l = {Ti-l)(l + T l 
TV 



m > 0, 



1 = (i-r i )(^r 1 + r r 2 + --- + ^" 1 )> ™<o, 1 = 1,2, 

it follows from Theorem [l] that we can write the equations (10 11 12 131 as standard 
conservation laws. Thus, the assumption that a generalized symmetry exist implies the 
existence of some conservation laws. 

If we assume that a second generalized symmetry exists, i.e. we can find a nontrival 
function G such that 



u o.o,t — 5o,o — G{Uh,o, Ufi-i,0) • ■ ■ , Wfl',0) u £, u ■ ■ ■ , % j,,), 



(14) 



where h, h' , k, k' are its orders, then we can state the following Theorem: 
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Theorem 2. Let the PAE m,i — F possess two generalized symmetries of orders (n, n! , k, k') 
and (n, n' , k, k'), Moo.t = goo and u m i = goo > o,nd let their orders satisfy one of the following 
conditions: 

Case 1 : n > 0, n = n + 1 Case 2: n < 0, ri = n 1 — 1 
Case 3: k > 0, fc = jfe + 1 Case 4: k' < 0, fc' = fc' - 1 

TTien m correspondence with each of the previous cases the PAE m t i — F admits a con- 
servation law 

(Ti - = (T 2 - l)g^ } , m = 1, 2, 3, 4, (15) 

where 

(1) , /. (2) , /«o,o (3) i j- (4) i fv-o,o Mr\ 

Po,o = 1 og/« 1 . , Po,o = lo S-^— » 9o,o =log/»o,i' 9o,o =1 °g7 — • ( 16 ) 

So the assumption that the PAE u\^\ — F have two generalized symmetries implies that 
we must have four necessary conditions of integrability, i.e. ther e mu st exist some functions 



of finite range Oq q , Oq q , p q , p a q satisfying the conservation laws ( 15 ) with p Q 5 , p 5 , g q , q q 

defined by eq. Q. g and Oq 2 q may depend only on the variables m^, and p 3 Q and j>q 4 q 
on u 0j . 

Summarizing the results up to now obtained we can say that a nonlinear partial difference 
equation will be considered to be integrable if it has a generalized symmetry of finite order, 
i.e. depending on a finite number of fields. This provide some conditions which imply the 
existence of functions p "^ or of finite range whose existence is proved by solving a 
total difference. 

For a DAE, when all shifted variables arc independent the proof that a total difference 
has a solution depending on a finite number of fields, i.e. is a finite range function, is carried 
out by applying the discrete analogue of the variational derivative, i.e. a function q n is (up 
to a constant) a total difference of a function of finite range iff 

^ = £^%^ = o, (17) 

ou n t-r* ou n+j 

see, e.g. [3"T] . For PAE this is no more valid as the shifted variables are not independent 
as they are related by the nonlinear PAE, in our case Ui t i — F(uifi, Uo.Oj u o.i)- This turns 
out to be the main problem for the application of the formal symmetry approach to PAE. 

To get a definite result we limit our considerations to five points generalized symmetries, 
i.e. when : 

uo,o = .9o,o = G(u lfi ,u-ifl, uq.o, uo,i> m o,-i), fl , «i, D fl , «-i,oS'«o,iS'«o,-i 7^0. (18) 

The existence of a 5 points generalized symmetry will be taken by us as an integrability 
criterion. This may be a severe restriction as there might be integrable equations with 
symmetries depending on more lattice points. However just in this case we are able to 
get sufficiently easily a definite result and, as will be shown in the next Section, we can 
even solve a classification problem. In this case we can state the following Theorem, which 
specifies the results obtained so far to the case of five point symmetries: 

Theorem 3. If the PAE ui t x — F possesses a 5 points generalized symmetry, then the 
functions 

1o"o = Q (m) ( M 2,o,wi,o,uo,o), m = l,2, 
Pofl = ^ (m) ( u o : 2,w ,i,uo,o), m = 3,4, 



must satisfy the conditions {15, 16) 



Then, using the relations (10 13) with n = k = 1 and n' = k! = —1, we get the fol- 
lowing relations between the solutions of the total difference conditions and the generalized 
symmetry G: 

<lo,o = ~ T i lo g G ,«i,o; 1o 2 o = T i logG„_ 10 , 

Po 3 o = ~ r 2 log G uo , ptfl = T 2 log G jM0 _ x . 
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So, to prove the integrability, which for us means find a generalized 5 point symmetry, for 



a nonlinear PAE u\\ = F, we have to check the integrability conditions (15 16). If they 



are satisfied, i.e. there exist some finite range functions q^ and Poo , we can construct 
the partial derivatives of G. The compatibility of these partial derivatives of G, given by 



eqs. (20 1, provides the additional integrability condition 



G 



G 



"0,-1, "0,1" 



(21) 



If these additional integrability conditions are satisfied, we find go,o up to an arbitrary 
unknown function of the form f(uo,o); which may correspond to a Lie point symmetry. 
This function can be specified, using the determining equations ([8]). 

The 5 point generalized symmetry go,Oj so obtained, will be of the form: 



.9o,o = $(«i,0j "o,o, «-i,o) + *(«o,i, "o,o, "o,-i) + v(uo,o). 



(22) 



3. Application of the test: an example 



Qoo* ( TO = 1, 2) and Po™ J (m = 3,4). This is not an easy task even if they are linear first 
order difference equations. A solution always exists but nothing ensure us a priory that 
the solution is a finite range function. So let us present a scheme for solving explicitly 
the integrability conditions we found for the equations on the square i.e. for finding the 
functions q$, q$, p^ and 



To check the integrability conditions (15 16 ) we need to find the finite range functions 



As an example of this procedure let us consider the solution of eq. (15) for m = 1, where 
Pofi = lo g(/tn,o)j %fi = Q^\u 2< o,ui,o,u 0t o), T 2 q^l = Q (1) (u 2 ,i,ui,i,m ,i)- (23) 



In eq. (23) we have the dependent variables u 2 ,\ and u± t i where u 2 .\ = F(u 2 ^ 1 "1,0, Wi 1) 
while = F(uifi) u o,o, "o,i)- So eq. (15) for m — 1 will contain the unknow function 



F which characterize the class of equations we are considering twice, one time to calculate 
u\ t \ in terms of independent variables and then to calculate 1*2,1 in term of u\^\ and of the 
independent variables. This double dependence makes the calculations extremely difficult. 
To overcome this difficulty we take into account that we are considering autonomous equa- 



tions which are shift invariant. So we can substitute eq. (15) for m = 1 with the following 
equivalent independent equations 



(m) 

Po,o 



P ( -10 = Q (m) (wi,l, 1tQ,l, «-X,l) - Q (m) (■"!,(), Wo.OjU-l.o) 



M 
Po,-i 



(m) 
■P-1,-1 



(24) 
(25) 



Q (m) ( u i,o,"o,o,u-i,o) - Q (m) (ui,-i,wo,-i,M-i ,-1) 
where, to simplify the notation, we introduce in the following the functions 

"1.1 = / (1,1) ( u l,0. u 0,0,f0,l) 1 "-1,1 = / ( ~ 1 ' 1) ("-l,0,'«0,0>"0,l) ; 

"1,-1 = /^ 1 ' _1) ('"l,0,«0 1 0,W0,-l), U-1,-1 = /( _1 ' _1 )(u_i j o,li0,0,«0,-l)) 

to indicate /o,o and its analogues. Moreover, we introduce the following two differential 
operators 



A 



B 







«0,0 



J «o,o a 



f (-14) 



a,. 



f (-l,l) " u -i,o' 



(26) 



f (i-i) f 

JU0,0 n JU0,0 



(-1-1) 



f (i,-i) " 10 



_1 _1) u U-l,o ■ 



A— 

Ju-l 



in such a way that the functional equations (24 25) reduce to differential monomials [T]: 
^W( My)B0|liU _y) = 0, BQW(tti,_i,uo,_i,u_i,_i) = 0, (27) 

AQW(u 1<0 ,u ,o,u- li0 ) = r( m > 1 \ BO (m) («^.«o 1 o,«-i,o) = rM. (28) 



Eqs. (27) are, by construction, identically satisfied while eqs. (28) provide a set of equa 



tions for the derivatives of Q( m ) (ui.o, "0,0, u-1,0) with respect to its three arguments. By 



8 



D. LEVI AND R.I. YAMILOV 



commuting the two operators ( 26 1 we can obtain a third equation for the derivatives of 
o, uo,o, u -i,o) with respect to its three arguments: 



[AS]Q (m) (wi,o,wo,o,"-i,o) = r (m ' 3) . 



(29) 



Eqs. (28 29 ), if independent, define uniquely the derivatives of the function Q( m ) (1*1,0, u o.o, u -i.o) 
and, if their consistency is satisfied, from them we get the functions themselves. 

In a similar manner from (T\ — l)Po"o^ = (-^2 ~ -O^o with m = 3,4 we get the function 

Poo = (uo,2i u o,ii Uo t o) and consequently the symmetry (22). 

This procedure works if the function F is known, i.e. if we check a given equation for 
its integrability. It also works if F is known up to some unknown arbitrary constants to be 
specified. In such case we solve a classification problem with unknown constants. However, 
the problem is much more difficult if F depends on unknown arbitrary functions of one, 
two or three variables. In such a case the coefficients of the operators (26) and functions 
r( m ' fc ) will depend on unknown functions, and A" 1 '^ may even depend on the composition 
of unknown functions. In this case a more complicated procedure might be necessary. 

3.1. A concrete example. Let us consider the following PAE [TU] 

2(u ,o + + «i,o + «o,i + 7[ 4u o,o«i,i + 2tti,oUo,x + 3(u ,o + ui,i)(ui i0 + u ,i)] $$0) 

+ (£2 + U)^0fiUl,l(ui,0 + «o,i) + (6 - U)ux,ou 0l i(u fi + Ul,l) + 
+ C"o,o«i,i«i,o«o,i = 0. 



Eq. ( 30 ) is a dispersive multi-linear partial difference equation which passes the A3 multiple 

]. Applying the Mobious transformation Uij 



7) we 



scales integrability test 
can rewrite it in a simplified form as 

(wo.owi.i + a)(wi,o + u o,i) + (2wi,o%,i + P)(uo,o + + 5 = 0, (31) 

where a, j3 and 5 are well defined functions of 7, £2, £4 and £. We now apply to eq. 
(31) the procedure outlined at the beginning of this section. Eq. (31) depends on three 
free parameters and we look for conditions on the three parameters, if any, such that the 
equation admits generalized symmetries. We get that the conditions are satisfied only in 
two cases: 

(1) a = 2/3 0, 5 = and, as (3 7^ we can always set (3 — 1. This choice of the 
parameters a and f3 corresponds to £ 2 = 3^4 + 37 2 and C 
corresponding integrable PAE reads: 



(«o,omi,i + 2)(ui, + u ,i) + (2ui,otio,i + 1)(«q,o + = °- 
In correspondence with the eq. (|32|) we get the generalized symmetry 



I27C4 in eq. (|30j. The 

(32) 



U0,0;i 



«o-2)(2< -l) 



B 



M0,0 
1 



1 



1 



(33) 



1 



_"0,1 M 0,0 + 1 U -iU 0fi + l_ 

(2) [3 = 2a 0, S = and, as a ^ we can always set a = 1. This choice of the 
parameters a and f3 corresponds to £2 = 6"f 2 — 3£t and £ = 127(7 2 — £4) in eq. (30). 
The corresponding integrable PAE reads: 

(1 + u fiUi t i)(uifi + it , 1) + 2(1 + uq,iUi,o)(uq,o + = 0- 
In correspondence with the eq. ( [34] ) we get the generalized symmetry 

UQ.-l 



(34) 



U ,0;t = A(U Q , 



7 + B ( u o.o - 1) 



U-inUin - 1 



W0,-1«0,1 - 1 



(35) 



Here A, B are constant coefficients, and in both cases, (A = 0, B ^ 0) and (A ^ 0, B = 0), 
the nonlinear DAE (33 35) are, up to a point transformation, equations belonging to the 
classification of Volterra type equations done by Yamilov [5T]. This shows that the eqs. 
(32 34) do not belong to the ABS classification. Moreover this calculation shows that the 
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A 3 integrability in the multiple scale integrability test is not sufficient to select integrable 
PAE on the square having five points generalized symmetries. 
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